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Abstract 

The Kondo problem is studied using the unitary Lie algebra of spin-singlet fermion bilinears. 
In the limit when the number of values of the spin N goes to infinity the theory approaches a 
classical limit, which still requires a renormalization. We determine the ground state of this renor- 
malized theory. Then we construct a quantum theory around this classical limit, which amounts to 
recovering the case of finite N. 
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I. INTRODUCTION 



Renormalizability-especially asymptotic freedom- is a deep and fundamental property 
of four-dimensional non-abelian gauge theories[l]. Unfortunately, at the moment we under- 
stand these theories well only in perturbation theory. Although qualitative connections with 
AdS theories of gravity have been made[2j, a precise quantitative solution for the spectrum 
of a non-abelian gauge theory remains out of our grasp. If we could get a non-perturbative 
formulation of the theory, free of divergences, we would be closer to the goal of proving that 
there is a mass gap and string tension[3j. 

The Kondo problem |1] provides us with a simpler, but quite deep, example of 
renormalization-and asymptotic freedom- with its own intrinsic physical interest. In a sense, 
it is the "hydrogen atom" of renormalization theory. It should be possible to solve this basic 
problem by simple methods, without unwieldy numerical calculations or clever ansatzes that 
depend on too many details of the particular model. Once we understand how to do this, 
we can see how similar ideas might apply to gauge theories and their dual gravities. 

Electrons in a metal move more or less like free particles, occasionally scattered by the ions 
which oscillate around their equilibrium positions by thermal fluctuations. As the temper- 
ature decreases, the ions move less and the resistance should decrease. It does, except that 
at some low temperature (T ~ 10K) the resistance starts to increase again, rising to a finite 
value as T — > 0. Kondo's explanation was that metals can have magnetic impurities (e.g., 
Iron atoms embedded within a Copper lattice) whose magnetic moments become ordered at 
low temperatures. These little magnets can scatter electrons too, an additional contribution 
to resistance at low temperature. The size of an atomic impurity is very small compared to 
the wavelength of the electron: it is much like a 5 function interaction. This causes di- 
vergences: Kondo calculated the magnetic contribution to resistance to be proportional to 
logT -1 . Thus Kondo explained why resistance grows at low temperatures. Theflaw is that 
it predicts infinite resistance at zero temperature, contradicting experimental observations. 
This is the Kondo problem. 

Wilson's landmark solution[5] follows a numerical approach (the Numerical Renormaliza- 
tion Group). The essential complication is that there are an infinite number of electrons, all 
of whom contribute to the divergent part of the interaction. The ingenious methods Wilson 
devised have not been generalized to the case of multiple impurities. Also, they have not 
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yet been made mathematically rigorous. 

Andrei and Wiegmann solved [HI U\ the Kondo problem using Bethe Ansatz meth- 
ods. While very elegant mathematically (and potentially rigorous), the method is not flexible 
enough to be useful beyond the basic example of the single impurity and with a linear dis- 
persion relation : in essence it is a fiendishly clever guess that happens to work because of 
the special symmetries of the problem. 

Nozieres[8j introduced a simple physical picture that is quite appealing: the magnetic 
impurity forms a bound state with an electron (or hole) and then mostly decouple from 
the rest of the electrons, leading to a Fermi liquid. It would be great to recover this from 
the more fundamental Kondo hamiltonian: somewhat analogous to recovering the chiral 
model from QCD. Affleck and Ludwig [9jcome closest to recovering such a picture based 
on conformal field theory: representations of Virasoro and Kac-Moody Lie Algebras. The 
motivation of our current work is in part to extend their studies beyond the critical point, 
to situations where conformal symmetry may not be exact. 

We work with a larger Lie algebra than Kac-Moody or Virasoro, spanned by all the 
spin-singlet bilinears of the electrons: even those that are not local field operators. We will 
see that these observables have small fluctuations in the limit where the spins take a large 
number of values N: the bilinears have commutators of order j|. Thus, even as h is kept 
fixed (for example, set equal to one) the commutators can be approximated by classical 
Poisson brackets. This "neo-classical" limitjin] retains many of the essential features of the 
theory, like the logarithmic divergence of the coupling constant. We will be able to perform 
renormalizaton explicitly and obtain the ground state in this limit. We then quantize the 
excitations around this ground state, thus recovering the finite N theory. This strategy was 
applied earlier to two dimensional QCD [[12], (which is free of UV divergences) as well as 
the Chiral Gross-Neveu model[13j(also called the non-abelian Thirring model). The Kondo 
problem is in fact simpler than either of these cases. Nevertheless, we believe it is instructive 
to work it out explicitly. 

The large N limit of the Kondo problem has been studied before [IT]. But the usual 
diagrammatic approaches are not the best way to understand how a new non-perturbative 
ground state forms. Our methods are more similar to the variational principles of BCS 
theory: but they allow for a systematic expansion around the neo-classical answer. 

While there remain technical and mathematical details to be explained, we recover a 
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simple picture: an impurity-electron condensate forms, and the excitations around this 
ground state are electron-like quasi-particles, but with a modified spectrum of energies. 

There has been a revival of interest in the Kondo effect because it occurs in quantum 
dots. Since the parameters can be tuned[14], this experimental realization holds promise 
of testing the theory as well as the potential for the invention of new devices based on the 
Kondo effect. 

II. FERMION BILINEARS 

Define fermionic operators satisfying the Canonical Anti-Commutation Relations, 

[A^,A 1e ,) + = 8f^ 
[A^A da ,} + = 6 a a , 

all other pairs of anti-commutators being zero. Throughout this paper (even where we 
speak of classical dynamics and Poisson brackets), we will use units such that H = 1; in 
particular, it is not equal to zero. Here, k labels the momentum of a conduction band 
electron, d an impurity state, and a the spin. Although a takes just two values in the real 
system, it will be convenient to let it take N values. There is only impurity, d — 1 . (We 
hope to generalize later to the case where several impurities are present.) Capital letters 
K, L will denote indices that can take either conduction band or impurity values: K = k 
or d. The momenta take a finite range of values k = —A, • • • , —1, 1, ■ • • A. The interesting 
physical region is when the energies are small compared to A, or equivalently, A — > oo. But 
this limit is very subtle, requiring a renormalization of a coupling constant. Our strategy 
will be to understand this first in the case of large N, and only then pass to the case of finite 
N. 

Define the spin-zero bilinears 

= ^A La . 

It is straightforward to check the commutation relations (the idea that fermion bilinears 
span a unitary Lie algebra go back to Schwinger in the early days of quantum field theory.) 
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= If [8?*5-s%*?]. 

Thus, in the limit of large N, these commutators become small: the quantum fluctuations 
in the spin-zero bilinears are small. They will tend to classical observables, the commutators 
being replaced by Poisson Brackets. We can think of the finite N the quantization 

of these Poisson Brackets, with playing the role of H in the usual quantum theory. This 
approach to the large N limit is motivated by the theory of solitons in the theory of strong 
interactions. In the large N limit, the Heisenberg equations of motion for the bilinears 
will tend to Hamiltonian equations. Such unusual classical limits (where a classical theory 
emerges even as h is finite, but some other parameter goes to zero) have been called 'neo- 
classical' in another context. We will solve for the static solution of least energy. For the 
solution to remain well-defined as A — > oo, we will have to renormalize a coupling constant 
even in this neo-classical limit. 

Then we have the P.B. in the large N limit 

-* » *JV ) = \pL ® N ~ d N § L \ . 

III. THE KONDO HAMILTONIAN 

The Kondo hamiltonian is 

H = "kA^Ato + JA^A da A^ A k , a , 
k k k' 

The sum over all momentum of the conduction band electrons amounts to evaluating the 
operators at the position of the impurity (the origin). 

We have allowed the spin generators to form a U (N) Lie algebra rather than SU (N) : 
this is a minor change, as the added singlet decouples from the rest. But it makes for easier 
book keeping. (In 't Hooft's study of the large limit of gauge theories, a similar passage 
from SU(N) to U(N) gauge group is made.) 

Of interest is anti-ferromagnetic case J > 0, where the electron and the impurity will 
form a spin-singlet. Define the operators 
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which involve the mixing of a conduction band state and a state located at the impurity. 
Expressed in terms of the bilinears, the Kondo hamiltonian is 

H = J2 - J&M = H - JH, 

k 

Assume that the energies of the conduction band electrons are u k are non-degenerate: 

iu k if k^l 

and that 

LV_ k = —UJ k . 

(Charge conjugation symmetry.) Moreover, we assume that the dispersion relation is asymp- 
totically linear: 

hm — — c 

|fc|->oc k 

for some constant c. Unlike in the Bethe Ansatz method, our approach does not rely on an 
exactly linear dispersion relation. Note that we do not allow u k ov k to take the value 0; this 
is to avoid an annoying zero mode and is not an essential restriction. 

IV. THE EFFECTIVE HAMILTONIAN 

Define the effective hamiltonian to be the matrix 



Then 



h\ = uj k 5? 

h k d = -9($), K = -g*m (1) 
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where, 



Using the fact that the Poisson brackets are those of the Unitary Lie algebra, it is easy 
to check that the equations of motion in the large N limit are, in matrix language, 

Note, by the way, that $^ is a conserved quantity. 

In particular, a static solution (such as the ground state ) will satisfy 

[/*(*),*] = 0. 

V. DIAGONALIZATION OF THE EFFECTIVE HAMILTONIAN 

For a given value of g, the eigenvalue equation for h is, 

(u k -u)U k -gU d = 



So that 



-g*U* = vU d 



\o\ 2 

(u k - v)U k + ^-U % = 



V 



and 



U k = \g\ 2 —^ -U*. 

v[y - LO k ) 

Summing over A;, the factor u' cancels out. The eigenvalues are then determined by the 
roots of the characteristic function 
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K 

Since the u)k are odd in k, this characteristic function is odd as well: 



X (-u) = -X{u). 

Thus there is always a root v — 0; the remaining roots appear as pairs differing by a sign. 
It is useful to use this symmetry to combine the k and —k terms in sum and write it as 



X{y) = v 



k=l < - u2 



For each root u a the eigenvector is given by 



U. 



uj k — v a 



This eigenvector will have length one if we set 



\U, 



d\2 



since 



X'(u) = 1 + E 



\9\ 



\2 • 



k {v-Lu k y 

Thus we have a 2A + 1 dimensional unitary matrix u% that diagonalizes the effective 
hamiltonian: 



h = f/diag(z/)f/ t . 



VI. THE GROUND STATE 



Since the static solution satisfies 



it must also be diagonalized by U : 



* = f/diag( / u)f/ t 
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The eigenvalues /x Q of are determined by the condition that it describe the ground state 
of the system: the negative energy states are occupied and the positive energy state is empty. 
The zero energy state carries any electrons that are left over after these assignments: 



1, 

0, 



v a < 



u a = 
v a >0 



The parameter fj, is the total number of electrons divided by N, modulo one. If /iq — \ 
we have just the right number of electrons to have a ground state that is invariant under 
the charge conjugation symmetry. It is useful to use instead a parameter that measures the 
departure from this symmetric case: 

i = \*-\ 

If N = 2, an odd number of electrons correspond to £ = and an even number of electrons 
to£=! 
Thus 



_ 1 - sgn(u a ) 

l^a — ~ r C,O a ,0- 



Of special interest are the elements 



*5 = X) U ^- 



Separating out the zero-mode contribution that is not charge conjugation invariant, the 
remaining sum can be written as a sum of residues: 



9t 



\9\ 



+ 9: 



dz 



2ni J D X(z) (u> k - z) 

Here D is a contour that starts at infinity a bit below the negative real axis, goes through the 
origin (where a principal value is taken) and then goes to infinity a bit above the negative 
real axis. The only poles arise from the roots of x( z ) a t which the residues are given by 
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the sum above, plus a principal value contribution from the origin. We can now deform this 
contour to go along the imaginary axis, and get the formula 



-f 



dz 



2m J D X(z) (uj k - z) 2ix J-ooX^y^ujk-iy) 



r 

j — c 



dy 



Thus, 



where 



9 r dy 



2ttJ X(iy) 



1 



= 9 -v[ 

k Jo 



uj k -iy uj k + iy 



dy iy 



o ,4 + ^X(iy) 



k 9^ 



u k 



E 



+ 



>-■? r 

k Jo 



dy 



[cul + yi] [l + %|*E(i/)] 



m>0 



This sum converges even as A — > oo; also the E(y) ~ ^ f° r large y. 
Setting 



V = \Wk\x, 



we get for the part even in k, 



,k 9 



f 

Jo 



dx 



® +d 7i\uj k \J [l + z 2 ][l + 2|<?| 2 E(K|:r)] 



It follows that 



This will be useful for renormalization. 

It is possible to get evaluate the sums when the spectrum is exactly (not just asymptot- 
ically) linear u k = ck . 



E(y) = 



-c + vrycoth 

2q^ 
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Figure 1: -3 -2 -1 1 2 3 

The impurity-electron condensate when £ = 0. Note the symmetry k — ► —k. 



2.0 r 



1.5 



1.0 



0.5 



Figure 2: -3 -2 -1 1 2 3 

The impurity-electron condensate when £ = 0.5. The symmetry k — > —k is broken. 



With 



we have 



a 2 da 



[1 + C(2)|^'| 2 ] fc ir J [k 2 + a 2 } [a 2 + \g>\ 2 (7racoth(Tra) -1)} 
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VII. RENORMALIZATION 



So far we have studied the problem for a fixed value of A. If we look back at the expressions 
in the last sections, we see that the limit A — > oo is convergent for S, A, \l> and h, as long 
as g is kept fixed. Since 



= JJ>? 



k 

and the sum is log divergent, it follows that J ~ T^K- More precisely, (remembering that 
only the even part of ^ contributes to the sum over k) 



A 1 

lim J (A) V — = 1. 

k=l 

This is asymptotic freedom. During renormalization we trade the divergent constant J~ l 
for g , which remains finite as A — > oo. In detail, 

j~u K ]= ^ir ^ 

{ ,9) foWo [<4 + y*][l + 2\g\*E(y)Y 
The dependence on g is sub-leading order in A. 

It is worth noting that g is a complex-valued parameter, although the original anti- 



ferromagnetic coupling J is real. The symmetry generated by the conserved quantity $ 



(the number of electrons occupying the impurity site) is spontaneously broken, as it corre- 
sponds to the phase of g. If there is a lattice of impurities, this would become a translation 
invariant field that breaks the gauge invariance of electromagnetism spontaneously: a pos- 
sible mechanism for superconductivity in heavy fermion systems [16j. 

VIII. THE RENORMALIZED THEORY 

The sum defining the characteristic function is convergent in the limit A — » oo : 

2M 2 



Removing the overall factor of i/(which just gives the obvious root at v = 0), we get the 
function 
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The characteristic function xi(^) f° r a nearly dispersion relation. 



Xl (i/) = i+£; 



2M 



In each interval [u>k,ujk+i} (with cjfe > 0) Xi( u ) increases from — oo to oo monotonically. 
Thus it has exactly one root u a in each such interval. 

The ground state \l/ and the effective hamiltonian h = hfif) of Q also make sense as 
A — > ookeeping g fixed. 

We can use the departure from the ground state as the dynamical variable : 

$ = _|_ 0. 

In addition, it is natural to rotate to the basis in which h(SSf) is diagonal. That is, put 

<p K L = 

and use the components 0^as our dynamical variables. Then the effective hamiltonian 
becomes 



= u a fi$ - [r(<f>)U*Uf + h.c] 

where 
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r(0) = jJ2€u*/u;. 

7,5 

Now, 



11* = -v — 

9 



-j y^^a' 



n * n * / j 



9* 9* y/X'(u a )X'M 

Because of the overall factor of J, this will vanish unless the sum in a, (3 diverges. Thus, 
for those 0^with just a finite number of non-zero entries 

W) = vJ% 

Such finite rank configurations satisfy the linear evolution equation equation 

d( t> a a r 



dt 

This subset is closed under time evolution. The corresponding quantum states, are free 
quasi-particles. 

IX. NEW DYNAMICAL VARIABLES 

The dynamical variables 4>p obtained after subtracting the static solution \l>, and passing 
to the basis diagonalizing h(ty), satisfy the Poisson brackets 

- i {</%, 01} = 5}^ - 8W} + - M 7 ) 8f% (2) 

(Recall that /i a are the eigenvalues of the static solution \l>.) This is the central extension 
of the unitary Lie algebra, defined for example, in the book by Pressley-SegalpJa]. If only 
a finite number of the 0^ are non-zero, we can supplement this with an element describing 
time evolution 

-i {h,<f>^} = [e a - e^. 
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All the effects of the impurity-electron interaction are contained in the shift of the energies 
from uj k to u a and in the occupation numbers \i a . The unitary transformation U% relates 
the new degrees of freedom to the old. In the continuum limit, this can be expressed as a 
scattering phase shift of the electrons. 

X. THE CASE OF FINITE N 

Now we are ready to return to the case of finite N. Since plays a role analogous to 
U, this amounts to quantizing the Poisson brackets (j2j. That is, find operators that satisfy 
these commutation relations 



\ - W + \Ma ~ fh\ ( 3 ) 



N 

The representation of interest is 

$ = ^ : a^ap, : . 

where a, are fermionic operators and the normal ordering is with respect to the Dirac 
vacuum of the energies e a : 

o) a | 0) = 0, e a < 

a a | 0) = 0, e a > 0. 
The hamiltonian just describes quasi-particles with these energies: 

a 

We get free particles only because we ignored terms in the hamiltonian that are not 
divergent. If we add UV finite interactions to the hamiltonian in addition, we get a Fermi 
liquid. 
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